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Abstract 

We report data on the thermal conductivity k(T, H) along the basal plane of the hexagonal 
crystal structure of superconducting Mgi_ y Alj,B2 with y = 0.02 and 0.07 at temperatures between 
0.5 and 50 K and in external magnetic fields between and 70 kOe. The substitution of Al for Mg 
leads to a substantial reduction of the heat transport via electronic quasiparticles. The analysis of 
the k(T, H) data implies that the Al impurities provoke an enhancement of the intraband scattering 
rate, almost equal in magnitude for both the a- and the tt band of electronic excitations. This is 
in contrast with conclusions drawn from analogous data sets for material in which carbon replaces 
boron and where mainly the intraband scattering rate of the a band is enhanced. Our complete 
data set, including new results of measurements of the low-temperature thermal conductivity of 
pure MgB2, confirms the validity of the Wiedemann-Franz law for both pure and doped MgB2- 
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I. INTRODUCTION 



In multiband superconductors, the onset of superconductivity may introduce gaps of 
different size on different parts of the Fermi surface. As a consequence, the thermal con- 
ductivity of multiband superconductors at T < T c , the critical temperature, and in external 
magnetic fields H < H C 2, the upper critical field, reveals specific features that are usually 
not observed for common type II superconductors. These features were first reported for 
MgB 2 (ref. []]) and subsequently also identified for NbSe2 (ref. [2]) and PrOs4Sbi2 (ref. []). 
The thermal conductivity k(H) of these materials in the mixed state of superconductivity, 
measured at constant temperature, exhibits an unusual two-step increase with increasing 
magnetic field between if cl and H c2 . The two steps reflect the separate recovery of the heat 
transport carried by quasiparticles of two electronic energy bands with significantly different 
energy gaps in the superconducting state.- 1 ^ In order to observe this type of behaviour of 
K e (H), the electronic contribution to the total thermal conductivity, it is important that 
in the corresponding temperature and magnetic-field ranges, the heat transport carried by 
phonons and represented by n ph (H), is relatively weak. In addition, also the conditions 
of weak interband scattering of the quasiparticles and a scattering of quasiparticles in the 
band developing the smaller gap compatible with a moderately dirty limit situation, need 
to be fulfilled. 4 All these conditions are ideally met for MgB 2 at helium temperatures and 
the two step enhancement of n e (H) with increasing field is very pronounced, such that the 
contribution of each of the two types of band to K e can be identified. Hence this allows, in 
principle, to investigate the influence of different perturbations on the heat transport car- 
ried by different energy carriers separately. Indeed, results of measurements of the thermal 
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conductivity of Mg(Bi_ x C x )2 (Ref. |6|) revealed that low concentrations of C replacing B 
reduce the step-like enhancement of k(H) near H c2 much more than the analogous feature 
of k(H) which is observed at much lower fields. This result suggests that the carbon related 
defects scatter quasiparticles of the band with the large gap much more than quasiparticles 
excited in the band with the smaller gap. This observation agrees with the expectation that 
C substitutions for B mostly enhance the scattering between states in the a band which is 
predominantly formed by boron sp x p y orbitals. In this work we describe and discuss results 
of analogous measurements of k(H,T) of Mg 1 _, / Al 2/ B 2 . The investigations of single crys- 
talline material with y = 0.02 and 0.07 reveal that Al replacing Mg at the few at-% level 
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reduces the heat transport of quasiparticles in both the a- and the tt band in a more or less 
equal manner. 

Aluminium replacing Mg in MgB 2 is expected to influence the superconductivity of the 
pure binary compound mainly via two effects.- 1 ^ First, Al doping causes a reduction of 
hole-type carriers, thus reducing the corresponding density of electronic states (DOS) and, 
likewise, the size of both gaps A CT and A„-, thereby accounting for the reduction of T c . The 
second consequence is caused by out of plane distortions created by Al occupying Mg sites 
which enhance the interband scattering thus leading to a relative increase of the gap in the 7r 
band. Indeed, results of various experimental studies confirm this view and suggest that for 
y < 0.1, the replacement of Mg by Al leads to a gradual reduction of A CT without substantial 
changes in the value of A n M-^^^ For y > 0.1 the situation is more controversial but 
most experimental results are compatible with the reduction of both gaps upon increasing 
y2&22il2. Because choosing y > 0.1 often leads to inhomogeneities in the Al distribution, such 
as the segregation of samples into Al-rich and Al-poor regions,— 1 ^ we restricted our study 
to material with y < 0.1, for which the influence of Al is expected to be as described above. 



II. EXPERIMENT 



The preparation and characterization o 
in this study is described in detail in ref. 



the Mgi_j / Al 2/ B 2 single crystals that were used 
I4I The synthesis under high pressure allowed 
the production of single-phase, single- crystalline material up to Al concentrations y ~ 0.1. 
For our measurements we selected two samples, below denoted as A12 and A17, with Al 
contents, as determined by EDX, of y = 0.02 and 0.07, respectively. Since we also intended 
to confirm or refute our earlier conclusion of the non validity of the Wiedemann-Franz law 
for MgB 2 at helium temperatures (ref. Ill), we repeated the measurement of the thermal 
conductivity of binary MgB 2 , using two newly synthesized single crystals, denoted as MB1 
and MB2 below. All samples were bars with the approximate shape of a prism and typical 
dimensions of 0.9 x 0.3 x 0.06 mm 3 . The thermal conductivity was measured using the same 
standard uniaxial heat-flow method that we applied in our previous experiments on pure and 
C-doped MgB 2 .- In this set-up, one end of the sample is attached to a copper body serving 
as a heat sink and at the opposite end, a miniature Ru0 2 resistor, serving as the heater, 
is attached. The temperature gradient along the sample was monitored employing either 
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RuC>2 resistor thermometers below 4 K or chromel/Au+0.07%Fe thermocouples above 4 K. 
These temperature sensors were each attached to a 25 pm diameter gold wire and the wires 
were glued to the crystals with silver-filled epoxy, approximately 300 to 400 /im apart along 
the direction of the heat flow. The temperature difference along the sample was chosen to 
be of the order of 3% of the average sample temperature. For the subsequent measurements 
of the dc electrical resistivity p, a four contact configuration was used. In order to minimize 
the uncertainty of the geometrical factor in the comparison of the electrical and thermal 
conductivities, the described thermometer contacts were also used as voltage contacts for 
the resistivity measurements. Nevertheless, the possible systematic errors in the separate 
evaluation of k or p values caused by the uncertainty of the sample geometry were as large 
as ±20%, whereas the random errors for k and p were only ±1% and ±0.1%, respectively. 



III. RESULTS 



The electrical resistivity was measured between 2 and 300 K in zero magnetic field and 
H = 50 kOe, oriented along the crystallographic c direction. For both new specimens of 
pure MgB2, MB1 and MB2, p(T) is very similar to the data that we obtained in earlier work 
(ref. [ll) for sample MB0. Likewise the values of T c , the transition widths and the residual 
resistance ratios turned out to adopt similar values as were established in our previous work. 
Note, however, that the absolute resistivities of the new samples just above T c are somewhat 
lower than before (see Tabled]). In the temperature range between T c and approximately 
150 K, p(T) may very well be represented by p(T) = p^ + aT 3 . The most obvious variation of 
p upon replacing small amounts of Mg by Al is an expected increase of the residual resistivity 
Pq, already reached at 40 K, which is caused by an enhanced scattering rate due to the Al 
defects in the Mg sublattice. The temperature dependent component due to the scattering 
of electrons by phonons changes only very little. 

In Fig. [1] we display the thermal conductivity of Mg 1 _ J/ Alj / B 2 (y = 0.02 and 0.07) as 
a function of temperature in zero magnetic field and for H = 50 kOe, with H oriented 
parallel to the c axis. In zero magnetic field, the absolute difference in magnitude of k is 
more pronounced at elevated temperatures but with decreasing temperature to much below 
T c this difference is reduced significantly and is only marginal at temperatures below 3 K. 
Below 3 K, the temperature dependence of k is fairly well represented by k(T) oc T 3 . In the 
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TABLE I: Samples of Mgi_ y Alj,B2 studied in this work. Sample MBO, listed for comparison, was 
studied earlier (see ref . U) . 

y T c (K) p(40 K) p(300 K)/p(40 K) 
from p(T) cm) 



MBO 38.1 ±0.1 2.05 6.80 

MB1 38.3 ±0.1 0.71 7.49 

MB2 38.2 ±0.1 1.20 7.11 

A12 0.017 35.9 ±0.15 2.36 3.90 

A17 0.073 33.0 ±0.2 4.13 2.70 



field induced normal state, the relative difference in the k values persists to temperatures 
below 1 K where k(T) approaches a variation linear in temperature. 

In Fig. [2] we display typical curves for k(H), obtained for the two investigated samples at 
various constant temperatures. Below H c \ and above H C 2, the thermal conductivity varies 
only weakly with field. However, in the mixed state, i.e., between H c \ and H C 2, k(H) exhibits 
a substantial field dependence, as expected. While H c i cannot reliably be determined from 
our data sets because the k(H) curves are irreversible in the vicinity of H&, the upper 
critical field can be identified quite accurately from the abrupt change of slope of k(H), as 
is indicated by the arrow in the upper left panel of Fig. [3j The resulting curves of H c2 (T) 
are shown in Fig. [3j We emphasize, as in our previous work, that H c2 (T) obtained in this 
manner is a reliable measure for the temperature dependence of the upper critical field of 
the bulk of the sample, insensitive to surface superconductivity or minority-phase inclusions 
(ref. 111). The observed trend of the alteration of H C 2(T), its values and anisotropics (see 
inset of Fig. [3]) upon Al concentration, is consistent with the results presented in ref. 
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FIG. 1: Thermal conductivity vs temperature in the ab plane of single-crystalline Mgi_ 2/ Alj / B2 
(y = 0.02 and 0.07) in zero magnetic field (open symbols) and H || c = 50 kOe (solid symbols). 
The arrows indicate the critical temperatures in zero magnetic field. 

IV. ANALYSIS AND DISCUSSION 

A. Multiband feature of thermal transport 

In MgB2 the total heat transport is carried by itinerant lattice and electronic quasiparticle 
excitations and k is given by the sum 
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FIG. 2: Thermal conductivity in the basal plane of Mgi_j / Al 2/ B2 (y = 0.02,0.07) vs at several 
fixed temperatures. The closed and open symbols correspond to the field directions perpendicular 
and parallel to the c axis, respectively. The arrow in the upper left panel demonstrates how H C 2 
for a given field orientation and at a given temperature was established. 



with n e = K e ^+K e ^-, where the last two terms capture the heat transport by quasiparticles of 
the a- and the n band, respectively. In principle, fitting the zero-field k(T) data in the tem- 
perature range T < T c by invoking eq. ([T]) and the appropriate theoretical expressions for the 
individual terms, estimates of the corresponding gaps in the quasiparticle excitation spectra 
can be made. The theory of Bardeen, Rickayzen and Tewordt (ref. Il6l ) has successfully been 
used to evaluate the energy gaps of conventional, single-gap superconductors from thermal 
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FIG. 3: The upper critical fields H C 2(T/T C ). The closed and open symbols correspond to the field 
directions perpendicular and parallel to the c axis, respectively. For comparison, our earlier results 
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for pure MgB2 (Ref. Il5l ) are shown by the dashed lines. The inset shows the anisotropies of the 



upper critical fields 7 = H^ 2 /H^ 2 - 



conductivity data. For multi-gap superconductors, more than one channel of quasiparticle 
heat transport in different energy bands need to be considered (ref. [171 ) and therefore, the 
number of fit parameters increases according to the increasing number of bands involved. If 
the heat transport carried by phonons is not negligible, the number of independent parame- 
ters increases even more. In particular, also the scattering of phonons by quasiparticles has 
to be taken into account. In our case, the upturn of k(T) in zero magnetic field just below T c 
indicates that the phonon contribution to the thermal conductivity in the normal state, K™ h , 
is not small in comparison with the normal-state component k™. Thus K p h is not negligible 
at T < T c and any fitting procedure in the manner indicated above seems unreasonable. 



Nevertheless, in ref. 
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an attempt of such an analysis has been made for data obtained for 



polycrystalline MgB2 where the phonon contribution to the total thermal conductivity is 
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claimed to be negligibly small. 

Information about details of quasiparticle scattering can be obtained from the analysis of 
the magnetic field dependence k(H) of the thermal conductivity, measured in the mixed state 
of type II superconductors at constant temperature. Usually, upon increasing H above H c i, 
k(H) first decreases because of an additional scattering of both phonons and quasiparticles 
by vortices. With still increasing H, the increasing density of quasiparticles in normal 
vortex cores enhances K e . Also for MgB2 the thermal conductivity in the mixed state is 
dictated by such processes but the resulting features of k(H) are clearly different from 
those observed for common superconductors. First, K e increases very rapidly in rather low 
fields, such that a considerable part of the normal-state thermal conductivity is recovered 
already when H <C H c2 . This is distinctly different from observations in common type II 
superconductors where, depending on the sample purity, K e (H) either grows significantly 
and rapidly only if H approaches H c2 (pure limit) or reaches its normal state value in 
an approximately linear fashion, i.e., K e cx H/H C 2 (dirty limit).— Second, when measured 
along different crystallographic directions and taking into account the anisotropy of the 
upper critical field, the low-field increase of k{H) does not scale as H/H C 2 as is the case 
for common type II superconductors. For MgB 2 , these anomalous features find a natural 
explanation by considering the anisotropy of the gap in the quasiparticle excitation spectrum, 
i.e., the different gap values for the a- and the n band, respectively. It turns out that the 
smaller gap is quenched already at H* < H C 2, where the critical field H* is practically 
independent of the orientation of H, reflecting the three-dimensional nature of the ir band. 
In this scenario, in the intermediate field regime H w H* < H C 2, the contribution of the 
quasiparticles created in the ir band is close to that in the normal state whereas in the a 
band, the large gap prevents the excitation of quasiparticles until H approaches H C 2 and 
hence, K e/J grows substantially only when H is close to H c2 . In shape this is similar to what is 
observed for common, moderately pure type II superconductors. Considering that at a given 
temperature, is significantly larger than H£ 2 , the ratio K e (H _L c ~ H*)/n e (H || c > H C 2) 
provides a rough estimate of the ratio k^^/k^. 

In Fig. H] we compare k(H), measured at T = 1.0 K, for pure and Al-doped MgB 2 ; 
only data points where k(H) is reversible are included in the plot. In this temperature 
range, K p h <C and therefore, the field dependence of k mostly reflects that of the thermal 
conductivity due to quasiparticles. The diagram clearly indicates that Al substitutions for 
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H (kOe) 

FIG. 4: Reversible thermal conductivity k normalized to its normal-state value K n vs. H at 1.0 K 
for pure and Al-doped MgB2- The closed and open symbols correspond to the field directions 
perpendicular and parallel to the c— axis, respectively. The inset shows doping dependence of n n 
at 1.0K. 

Mg reduce the ratio k^/k" and thus K™ n /K™ a but the latter ratio changes much less than 
the total electronic thermal conductivity in the normal state (see inset of Fig. Hj). For the 
ratios k^/k^ at 1 K, we estimate the values of 1.9, 1.6 and 1.2 for the samples MB1, A12 
and A17, respectively. This suggests that although the introduction of Al for B reduces the 
thermal conductivity of quasiparticles slightly more strongly in the ir band, the overall effect 
is nearly the same in both channels. This is in stark contrast to the observations made in the 
case of Mg(B 1 _ 2 .C x ) 2 where a 3% substitution of B by C reduces the quasiparticle thermal 
conduction of the a band to a negligible level in comparison to the 7r band contribution. 6 
The general equations capturing the thermal conductivity due to quasiparticles can be 
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simplified to the form 



(2) 



where Cy e is the electronic specific heat, vp the Fermi velocity, and l e the electronic mean 
free path. In the mixed state the mean free path is given by l e = 17^/(7 + 7^4), where 7 is the 
scattering rate due to impurities and 7^ is the scattering rate due to Andreev scattering of 
quasiparticles at vortices.— In this scenario, only 7 is affected by Al impurities. For MgB 2 
this approach has to be modified to 



with i = 71, a. It has been established that Al substitutions up to y — 0.1 change neither 
Cye nor vf,tt in the ah plane in any significant way— With an Al content y = 0.1, vf,o- in the 
ah plane is reduced only by 10 to 15% with respect to its value in the pure compound.— 1 ^ 
This suggests that the reduction of the quasiparticle mean free path l e is responsible for 
the substantial reduction of n e shown in the inset of Fig. HI If we assume that both the 
specific heat and the Fermi velocities remain indeed unchanged, our data imply that a 7.3 
at.% substitution of Mg by Al leads to reductions of and l e , n by factors of 2.7 and 
4.5, respectively. A de Haas-van Alphen study of Mgi_ ?/ Al 2/ B2 (ref. I20I ) indicates that for 
y ~ 0.075, the mean free path on the a sheet of the Fermi surface is approximately a factor 
of 2 shorter than for pure MgB 2 and for the tt sheet, the reduction of / ej7r is estimated to be 
of the order of 3 or slightly larger. These results are obviously in very good agreement with 
our conclusions. 

B. Relations between thermal and electrical conductivity 

The relaxation rates which govern the electronic charge- and heat currents in metals, 
respectively, are in general not identical. However, some of the involved scattering processes 
are, under certain conditions, equally efficient in relaxing both electric and thermal currents. 
In these cases, the electrical resistivity p and the electronic thermal conductivity K e are 
related via the law of Wiedemann and Franz (WFL), such that 



where L = 2.45 x 10 8 W Q K 2 is the Lorenz number. This relation is expected to 
be valid for any metal whose electronic subsystem can be described as a Fermi liquid. 




(3) 



Ke (T) = L T/p(T) 



(4) 
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Elastic scattering of electrons by static lattice defects is one of the mechanisms which are 
compatible with the validity of the WFL. If this type of scattering dominates all other 
processes, the residual resistivity at low temperatures adopts a constant value and, in the 
same temperature regime, K e varies linearly with temperature. The experimental Lorenz 
number L(T) = k(T)p(T)/T, with k(T) as the directly measured thermal conductivity, is a 
constant and equal to L$ only if the phonon thermal conductivity is negligible, i.e., K p h <C K e 
and hence n e « k. Inelastic scattering processes always lead to L(T) < Lq and sizable heat 
transport by phonons provides and excess heat conduction and hence L(T) > L . 

In the field induced normal state, all samples of MgB2 exhibit a temperature independent 
resistivity below 40 K, suggesting that defect scattering of quasiparticles dominates the low- 
temperature electronic transport. In passing we note that even in external magnetic fields, 
oriented along the c axis and of the order of 50 kOe, i.e., exceeding the corresponding bulk 
value of H c2 at T = K for all investigated materials, the onset of a superconducting 
transition is observed at very low temperatures. Most likely this transition is related to 
properties of the sample surface and hardly affects the predominantly probed heat transport 
in the bulk of the sample. In order to check the validity of the WFL, we thus used p = po ~ 
p(40 K), presented in Table fl] 

In part of our earlier work (ref. Q) we encountered evidence that for pure MgB2, 
L(T)/L > 1 at low temperatures which could not simply be attributed to a substantial 
phonon contribution to k(T), at least not below T = 6 K. Our most recent experiments on 
Mg(Bi_ x C x .)2 single crystals (ref. fj), however, gave no evidence for a violation of the WFL 
in the same temperature region. In Fig. [5] we plot these earlier results together with the 
corresponding new data that we accumulated for Al-doped MgB 2 . Below 1 K, the series of 
both C and Al doped materials exhibit ratios Lj Lq very close to 1. The smoothly increasing 
deviation from the ideal value with increasing temperature can easily be attributed to a non 
negligible K p h which, in the normal state of a metal exhibits a temperature variation that is 
typically stronger than linear in T. 

Since it appeared that the WFL is well satisfied in MgB2 with small amounts of defects 
due to Mg and B substitutions in the form of Al and C, respectively, it seemed of interest 
to verify whether the reported violation of the WFL is limited to pure MgB 2 . We thus 
reexamined the electrical and the thermal transport properties of two new single crystals 
of MgB2 (MB1 and MB2) between 0.5 and 6 K, both in the superconducting state in zero 
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FIG. 5: Normalized Lorenz number L(T)/Lq for three samples of pure MgB2 and for C- and 
Al-doped MgB 2 . C6 denotes Mg(B . 9 4C .o6)2 (data from Ref.M). 



magnetic field and in the normal state, induced by a magnetic field of 50 kOe, oriented along 
the c direction of the crystal lattice. The corresponding data resulted in the L/Lq ratios 
that are also displayed in Fig. [5l 

In contrast to the substantial and temperature-dependent enhancement of the Lorenz 
number that resulted from the data that were earlier obtained for sample MBO, we now 
note an almost temperature-independent L(T), somewhat smaller than L , for both new 
samples of pure MgB 2 . The reduction of L(T) to below the value of L is most likely not 
caused by inelastic scattering of the quasiparticles, because an almost parallel shift of L(T) 
from Lq would require that the responsible scattering would add a temperature-independent 
contribution to the thermal resistivity. However, all inelastic scattering processes that we 
are aware of, vary their strength with temperature. A more plausible explanation is to 
attribute the constant shift of L(T) from L to uncertainties in the sample geometry. If, 
as in our experiments, the same contacts are used to measure either the voltage or the 
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temperature difference along the sample in the respective measurements of the electrical or 
the thermal conductivity, the same geometrical factor applies for the corresponding equations 
for calculating p or k. Thus in calculating L(T), this factor drops out. This is really only true 
for cases, where the two contacts can be mounted far apart, in the sense that their separation 
d is much larger than the contact width S. Our samples are quite small and hence d ~ 5. In 
this case it is quite possible that the effective distances between the contacts for measuring 
the voltage or the temperature difference, i.e., d p or d K , are different. If now the ratio d K /d p 
differs from 1, the calculation of L will result in a number that differs from L . For our 
small crystals, the distance between the contacts was of the order of 300 to 400 pm and 
the contact widths varied between 50 and 70 /zm. It is therefore not unlikely that d K ^ d p 
which, in turn, may explain the constant shift of L(T) from L . 

In this sense we believe that no anomalous violation of the WFL is observed for the two 
new samples of single crystalline MgB 2 . This suggests that our earlier observation reported 
in ref. [l is not an intrinsic property of the material. We note, however, that the k(H) data 
exhibit the same features for all the investigated samples of pure MgBoand therefore, all 
the conclusions concerning the two-band features of MgB2 drawn in ref. Q remain valid. 

V. CONCLUSIONS 

Based on experimental results of the electrical and thermal conductivity we note that the 
substitution of magnesium by aluminium in MgB 2 provokes an increase of the quasiparticle 
scattering of similar magnitude in both the a- and the n band of electronic states. The 
calculated reductions of the quasiparticle mean free paths in the a- and 7r band, respectively, 
are in agreement with results of de Haas- van Alphen experiments. The considerable violation 
of the Wiedemann-Franz law at low temperatures for pure MgB 2 that we reported in ref. Q 



s of this compound, 
and this work), no 



could not be confirmed in additional experiments using new single crystal; 
Since also for all the doped MgB 2 crystals that we studied since (ref. 
anomalous deviations of the WFL were identified, our previous observation must have been 
of extrinsic origin. 
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